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In this paper we show quantum fluctuation effect of fully frustrated Ising spin 
systems. Quantum annealing has been expected to be an efficient method to 
find ground state of optimization problems. However it is not clear when to 
use the quantum annealing. In order to clarify when the quantum annealing 
works well, we have to study microscopic properties of quantum annealing. 
In fully frustrated Ising spin systems, there are macroscopically degenerated 
ground states. When we apply quantum annealing to fully frustrated systems, 
we cannot obtain each ground state with the same probability. This nature is 
consistent with "order by disorder" which is well-known mechanism in frus- 
trated systems. 
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1. Introduction 

Quantum information processing has been expected to be able to solve 
difficult problems which are not tractable on classical computer. ^"■^ Then a 
number of researchers have studied how to realize a quantum computer and 
quantum algorithm itself. Study on quantum information from a viewpoint 
of statistical physics came into the world about a decade ago. This is called 
quantum annealing, '^^-'^'^ in other words, quantum adiabatic evolution. In 
this paper we review on microscopic properties of quantum annealing with 
providing a specific example of Ising spin systems. 

The Ising model has been regarded as a standard model in statistical 
physics since it can be represented nature of phase transition and a couple 
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of physical phenomena in real magnetic materials. The Hamiltonian of the 
Ising model with random bonds is given as 

^ = - E -^^^^^i' < = (1) 

where (i, j) represents pairs of sites on the given graph. Since the spin 

variable al takes ±1, this variable obviously can be regarded as a "bit". 
Then the Ising model has been adopted not only for magnetic materials but 
also for information science. For example, traveling salesman problems and 
neural network can be mapped onto the random bond Ising model. These 
problems are categorized into optimization problems. It is difficult to find 
the best solution of optimization problems, since the number of candidate of 
solutions increases exponentially as the number of elements increases. Then 
we have to use engineered method to search the best solution of optimiza- 
tion problems. One of the methods is Monte Carlo simulation. However if we 
use Monte Carlo simulation, we face on difficulty of search the best solution 
of optimization problem, since energy landscape of optimization problems 
is complicated in general. In order to avoid such a difficulty, a number of 
researchers have improved the standard Monte Carlo simulation. 

One of the important methods is simulated annealing which was first 
proposed by Kirkpatrick.-'^^'^^ Roughly speaking, the simulated annealing is 
an optimization technique by decreasing temperature i.e. thermal fluctua- 
tion (Fig. 1). About fifteen years later since the publication of Kirkpartick's 
pioneering work, an alternative method of simulated annealing was pro- 
posed by Kadowaki and Nishimori.'* In order to obtain the ground state, we 
decrease quantum field as substitute for temperature (Fig. 1). The method 
is called quantum annealing. The quantum annealing is expected to be an 
efficient method to obtain the better solution of optimization problems. Ac- 
tually, the performance of quantum annealing has been demonstrated by a 
number of researchers. ^^"^^ 

However the solution obtained by the quantum annealing is worse than 
that obtained by the simulated annealing for some cases. It is not clear 
when to use and when not to use the quantum annealing. Then we should 
study microscopic mechanism of quantum annealing in order to understand 
why the quantum annealing works well for many cases. In the next section, 
we will review how to implement standard Monte Carlo simulation, simu- 
lated annealing, and quantum annealing. In Section 3, we will consider a 
quantum fluctuation effect of frustrated Ising spin systems. Finally we will 
summarize the microscopic mechanism of quantum annealing and show fu- 
ture perspective. 



January 27, 2013 0:51 



WSPC - Proceedings Trim Size: 9in x 6in ORC arXiv 



3 



Simulated annealing 

Temperature (T) 

Fig. 1. Schematic picture of simulated annealing and quantum annealing. The star 
denotes the position of T = and F = which is the target position. 

2. Computation Method 

Many optimization problems can be represented by random Ising models as 
stated in the previous section. The ground state of the random Ising spin 
system corresponds to the best solution of the given optimization problem. 
There can be many strategy to find a ground state of these systems. We 
often use Monte Carlo method as one of the most general algorithm. In 
the first part of this section, we will review on Monte Carlo simulation and 
simulated annealing. Next we will introduce quantum annealing which is 
regarded as an alternative for the simulated annealing. 

Here we consider the Hamiltonian given by Eq. (1). Algorithm of Monte 
Carlo simulation can be summarized as follows: 

step 1 By using random number generator, we prepare a random state as 
the initial state. 

step 2 We select a site at random and calculate an effective field. The 
effective field is defined as 

/^S= E J^^^h (2) 

J G n . n . oii 

where i denotes the label of the selected site. The summation takes 
over the nearest neighbor (n.n.) of the i-th site, 
step 3 We flip the spin of the i-th site according to some kind of transition 
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probability. How to choose the transition probabihty will be given 
below. 

step 4 We repeat step 2 and step 3 for a long time. 

There can be a couple of definitions of transition probability in step 3. For 
example, heat bath method and Metropolis method are very famous and of- 
ten adopted. Transition probability in heat bath method and in Metropolis 
method are given as 

w'hbK -o-f) = — TTTl TZUTT^ (^) 

[e ^''^"off {crtKii > 0), 

Both two transition probabilities obey the detailed balance condition. Very 
recently, an efficient algorithm was proposed by Suwa and Todo.^^ Their 
algorithm is not based on the detailed balance condition but on balance 
condition only. It helps us to obtain in a shorter time than conventional 
method, although it depends on problem. 

At high temperature, state can be changed easily since the canonical 
distribution is almost flat. However it is generally difficult to change state 
at low temperature because of complicated energy landscape and relatively 
high energy barrier. Then we can avoid partly such a difficulty by decreasing 
the temperature gradually. This is called simulated annealing. It has been 
widely adopted for many cases. 

Next we review how to implement Monte Carlo simulation for random 
Ising systems with transverse field. The Hamiltonian is given as 

H = He + ^qj (5) 

:= - E -^i^^i^h ■■= -r E (6) 

(hi) i 

where erf and erf represent the x-component and z-component of the Pauli 
matrix at the z-th site, respectively. Here we adopt transverse field as a 
quantum term. For small systems, we can obtain equilibrium physical quan- 
tities since; wo can calculate the Boltzmann weight c^^^ by using exact 
diagonalization. However it is difficult to calculate the Boltzmann weight 
Q-PH fQp relatively large system because of limitation of memory. Then, 
in general, we have to adopt other strategy for obtaining the equilibrium 
properties of large-scale quantum systems. Nowadays the quantum Monte 
Carlo simulation is regarded as an efficient algorithm. In order to use Monte 
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Carlo method, wc have to calculate equilibrium probability c~^'^ as men- 
tioned above. Here we calculate partition function of the original quantum 
system H by applying the Trotter-Suzuki decomposition.^^'^^ The partition 
function is given by 

Z = Tie-^^ = Tre-^^^'+^-^ = ^ (E|e-^(«=+«-)|S) . (7) 
By using the following relation 



we obtain 



Z= (Si|e-^|S;)(Sl|e- 

{Si,fe,*-,i,=±l} 



x(E2|e-^|S^) (S^le-^IEs) 



x(S„|e-^|Sj (E;„|e-'^|Si), 



where expresses direct product such as 



(9) 
(10) 



where N represents the number of sites and k denotes the position along 
the Trotter axis. The first subscript of a is the position in the real space. 
Since He is a diagonal matrix, then we obtain 



(11) 



On the other hand, we can calculate (E^|e as follows: 



(E',|e-^|Sfe+i) 



■ sinh 



\ m 



g[i iogcoth(^ Ef=i *i.fe*U+i)] (12) 



The partition function of the original Hamiltonian H is expressed as 







^20rv 


Z = lim 


- sinh 




m— >oo 


2 ^ 





{ai,fc=±l} 
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Then the equivalent classical effective Hamiltonian He can be written as 

= - E E t^u-h - is logcoth ^ E E (13) 

fc=l ^ ^ / i=l fe=l 

Note that ffi^m+i = ai^i because of Eq. (9). We can obtain equilibrium 
properties of the original quantum system H by calculating that of the 
classical effective Hamiltonian He- We can perform the quantum annealing 
by decreasing transverse field F. 



3. Quantum Fluctuation Effect of Frustrated Ising Systems 

In this section we consider quantum fluctuation effect of frustrated Ising 
systems comparing with thermal fluctuation effect. Although there are a 
couple of types of frustrated systems such as geometric frustrated systems 
and random systems, we focus on regularly frustrated systems in this paper. 

Figure 2 shows all ground states of three spin systems. The dotted lines 
and the solid linos in Fig. 2 represent ferromagnetic interactions and an- 
tiferromagnctic interactions, respectively. The white and black circles in 
Fig. 2 indicate + spins and — spins, respectively. As shown in Fig. 2, there 
are six ground states for antiferromagnetic case. The crosses of the right 
panel in Fig. 2 depict energetically unfavorable interactions. In general, 
there are many degenerated ground states in frustrated systems because of 
such unfavorable interactions. 

One of regularly frustrated systems is antiferromagnetic Ising model 
on triangular lattice. Figure 3 takes three of many ground states of this 
system. Figure 4 also takes three of many ground states of antiferromagnetic 
Ising model on kagome lattice as an another example. In these models, as 
the number of spins N increases, the number of ground states increases 
exponentially. 

6--h 

A AAA 




Fig. 2. The white and black circles denote + spins and 

The ground states of the ferromagnetic triangle cluster, 
of the antiferromagnetic triangle cluster. 



— spins, respectively. (Left panel) 
. (Right panel) The ground states 
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Fig. 3. Ground states of antifcrromagnctic Ising system on triangular lattice. The dot- 
ted circles denote free spins where the effective field is zero. The left panel shows one of 
the maximum free spin states. 




Fig. 4. Ground states of antiferromagnetic Ising system on kagome lattice. The dotted 
circles denote free spins where the effective field is zero. The left and middle panel show 
maximum free spin states. 

Our motivation is as follows: Suppose we consider systems where there 
are many ground states. When we apply quantum annealing to such sys- 
tems, which states arc selected at the final time? What arc similarities 
and differences between simulated annealing and quantum annealing? In 
order to consider them, we study quantum fluctuation effect of regularly 
frustrated Ising spin systems. 

Before we study quantum fluctuation effect, we consider thermal fluc- 
tuation effect for comparison. We apply simulated annealing to regularly 
frustrated Ising systems. Suppose the cooling speed is set to be very slow 
from high temperature to T = 0+. Then we can obtain each ground state 
with the same probability. This is a reasonable nature from a viewpoint of 
the principle of statistical physics. If the internal energy in the i-th state 
and that in the j-th state are the same, probabilities of both states are the 
same. This is called the principle of equal weight. Then such nature appears 
despite details of model such as lattice structure. 
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Simulated annealing i i Quantum annealing 




free spins 



number of 
free spins 



Fig. 5. Schematic picture of probability distribution at the final time for simulated 
annealing (left panel) and quantum annealing (right panel) for adiabatic limit. 

On the other hand, situation is drastically changed when we apply quan- 
tum annealing to regularly frustrated Ising systems even if we decrease 
transverse field very slow. The dotted circles in Fig. 3 and Fig. 4 represent 
spins where the effective field is zero. The states with the highest appearance 
probability corresponds to the maximum free spin states. As the number of 
free spins decreases, the appearance probability decreases at least in large 
number of free spin region. This result is consistent with the Villain model 
case which is a typical example of regularly fully frustrated Ising systems 
as well as antiferromagnetic Ising model on triangular lattice and kagome 
lattice. 

Finally, Fig. 5 represents a schematic picture of probability distribu- 
tion at the final time for simulated annealing and quantum annealing for 
adiabatic limit. 

4. Conclusion 

In this paper, we review implementation method of Monte Carlo simulation 
of classical Ising models and quantum Ising models. We also show micro- 
scopic properties of quantum annealing for fully frustrated Ising spin sys- 
tems comparing with that of simulated annealing. When we apply quantum 
annealing to fully frustrated Ising systems, we cannot obtain each ground 
state with the same probability. In contrast, we can obtain each ground state 
with the same probability by simulated annealing because of the principle 
of statistical physics. In this paper we consider adiabatic limit of quantum 
annealing. We have to study microscopic properties of quantum annealing 
beyond adiabatic limit keeping practical situations in our mind. 

In this study, we focus on transverse field as a quantum fluctuation. It 
is not necessary that we restrict quantum fluctuation to transverse field. 
Moreover we can adopt a new type of fluctuation such as invisible flue- 
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tuation which is first proposed by Tamura et al.'^^'^^ In order to improve 
original simulated annealing, we have to study fluctuation effect of random 
systems in a general way. 

Quantum annealing, in other words, quantum adiabatic evolution was 
born in frontier between quantum information science and quantum statis- 
tical physics. Study on quantum annealing is quite primitive and will have 
a favorable influence on both fields. 

The author is grateful to Masaki Hirano, Kenichi Kurihara, Yoshiki 
Matsuda, Seiji Miyashita, Hiroshi Nakagawa, and Issei Sato for fruitful 
discussion. The author also thanks Ryo Tamura for critical reading. The 
author is partly supported by Grant-in-Aid for Young Scientists Start-up 
(21840021) from the JSPS, MEXT Grant-in-Aid for Scientific Research (B) 
(22340111), and the "Open Research Center" Project for Private Universi- 
ties: matching fund subsidy from MEXT. The computation in the present 
work was performed on computers at the Supercomputer Center, Institute 
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